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Abstract 

We consider a planar differential system x = P(x,y), y = Q(x,y), where 
P and Q are C 1 functions in some open set and '=■§[■ Let 7 be a 

periodic orbit of the system in U. Let f(x, y) : U C M 2 — > E be a C 1 function 
such that 

df df 

where k(x, y) is a C 1 function in U and 7 C {(x, y) \ f(x, y) = 0}. We assume 
that if p G U is such that f(p) = and V/(p) = 0, then p is a singular point. 

We prove that ^ (if + dt = $ k ^)) dt > where T > is the 

period of 7. As an application, we take profit from this equality to show the 
hyperbolicity of the known algebraic limit cycles of quadratic systems. 
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1 Introduction 



We consider a planar differential system 

x = P(x,y), y = Q(x,y), (1) 

where P and Q are C 1 functions in some open set hi C M 2 , and ' = 4. A singular 
point of system ((TJ) is a point p E hi such that -PQo) = <5(p) = 0. We assume that 
all the singular points of ((TJ) are isolated. 

Given a system (0), we can always consider its vector field representation F(x, y) = 
(P(x,y),Q(x,y)). 

We will denote by div(x, y) the divergence of system ([TJ), that is, div = dP/ dx + 
dQ/dy. 

We also need to consider the flow of system (£Q), which we denote by $t(p) 
and which represents the unique solution of system (£Q) passing through the point 
p G hi C M 2 . We notice that for each p £ hi there exists an e p > (which may be 
e p = +oo) such that t G (— e p , e p ) is the maximal symmetric interval of existence of 
the solution of (JTJ) passing through p. We have that ^(p) = (P($ t (p)), Q($ t (p))), 
for all p G W and i G (— e p , e p ), and $o(p) = V- Given p E hi, the function $(•,£>) : 
(— e p ,e p ) — > M. 2 , where &(t,p) := &t(p), defines a solution curve or orte of (JTJ) 
through the point p. 

A limit cycle of system (P) is an isolated periodic orbit. Let 7 be a limit cycle 
for system (JTJ). We say that 7 is stable if there exists a neighborhood W 7 C W of 7 
such that for all p G W 7 , we have lim^ +00 <i($ 4 (p), 7) = 0. As usual, the previous 
application d is the distance between sets in the Hausdorff sense. Analogously, we 
say that 7 is unstable if there exists a neighborhood W 7 C hi of 7 such that for all 
p G W 7 , we have lim^.oo d($ t (p), 7) = 0. 

There might be limit cycles which are neither stable nor unstable. Using the 
Jordan curve theorem, which states that any simple closed curve, as the limit cycle, 
7 separates any neighborhood W 7 of 7 into two disjoint sets having 7 as a boundary, 
we can consider W 7 as the disjoint union of hii U7 Uhi e , where hii and hi e are open sets 
situated, respectively, in the interior and exterior of 7. When for any p G hii we have 
limt^ +00 <i($i(p), 7) = whereas for any q G hi e we have lim^^oo 7) = 

(or, the other way round, for any p G hii we have lim^.oo d(^ t (p), / y) = whereas 
for any q G hi e we have lim 4 ^ +00 <i($ 4 (g),7) = 0), we say that 7 is semi-stable. 

Any limit cycle 7 of a system is either stable, unstable or semi-stable as it 
is stated in [T^j. For a detailed description of the classical known results on limit 
cycles see also jT^j. 

The following result, which is stated as a corollary in page 214 of [TSj, gives a 
formula to distinguish the stability of a limit cycle. 
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Theorem 1 Let j(t) be a periodic orbit of system (Q) of period T. Then, 7 is a 
stable limit cycle if 

-T 

div(7(t)) dt < 0, 







and it is an unstable limit cycle if 



T 

div(7(t)) dt > 0. 

It may be stable, unstable or semi-stable limit cycle or it may belong to a continuous 
band of cycles if this quantity is zero. 

A sketch of the proof of this theorem is given after the forthcoming Theorem El 

When the quantity j Q div (t)) dt is different from zero, we say that the limit 
cycle 7 is hyperbolic. 

Since we are considering differential systems ((TJ in the class of functions C 1 , we 
may have a limit cycle 7 belonging to a sequence of periodic orbits {7„ , n G N} 
with 7„ + i in the interior of j n , such that the sequence accumulates to a singular 
point, a periodic orbit or a graphic and such that every trajectory between 7 n and 
7„+i spirals towards 7 n or 7„ +1 as t — > ±00. This kind of phenomena does not exist 
for analytic systems. 

In this work, we give another quantity which equals to J Q div(7(i)) dt for a pe- 
riodic orbit 7 defined in an implicit way, as explained below. This is the main 
result of the article and it is stated in Theorem El in the following section. We can, 
therefore, distinguish the hyperbolicity of a limit cycle using two different quantities. 

Given a planar system (JTJ) (or equivalently its vector field representation F(x, y) = 
(P(x,y),Q(x,y))), we define an invariant curve as a curve given by f(x,y) = 0, 
where / : U C M 2 — > M is a C 1 function in the open set U, non locally null and such 
that there exists a C 1 function in U, denoted by k(x,y), for which 

df Of 

P(%,y)-Q^(x,y) + Qfay) q-{x,v) = K x ^v) f( x ^y)^ ( 2 ) 

for all (x, y) G U. The identity (J2J) can be rewritten by V/-F = kf. As usual, V/ de- 
notes the gradient vector related to f(x,y), that is, Vf(x,y) = (§f y), y)), 
F(x,y) is the previously defined vector (P(x, y) , Q(x, y)) , and ■ denotes the scalar 
product. We will denote by / or by ^ the function V/ • F once evaluated on a 
solution of system ((TJ). 

We will always assume that if p G U is such that f(p) = and V/(p) = 0, then 
p is a singular point of system This is a technical hypothesis which generalizes 
the notion of not having multiple factors for algebraic curves. For instance, if we 
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had written that the periodic orbit 7 was contained in f 2 (x,y) = 0, then we would 
have that V(/ 2 )(p) = for all p G 7, in contradiction with the hypothesis. 

We notice that, as a particular case, we may have a function f(x,y) given by a 
polynomial in M[x, y). In such a case, f(x,y) = is called an invariant algebraic 
curve. When, in addition, the system is polynomial, that is, P,Q<E M[x, y), then the 
function k(x, y) is a real polynomial called cofactor. When we consider an algebraic 
curve, we can always assume that it is defined by a polynomial f(x, y) — such that 
the decomposition of f(x, y) has no multiple factors. The same assumption must be 
done for curves defined by C 1 functions and it is equivalent to the assumption that 
if p G U is such that f(p) = and V/(p) = 0, then p is a singular point of system 
P|). More explicitly, assume that the system (Q) has an invariant algebraic curve 
given by f(x, y) = and assume that the decomposition of f(x, y) in the ring WL[x, y] 
has no multiple factors, that is, f(x,y) = b\(x,y)b2(x,y) . . .bk(x,y) where bj(x,y) 
is an irreducible polynomial in M[x, y] and bi(x,y) 7^ cbj(x,y) for any c G K. — {0} 
if i 7^ j. Let p G U be such that /(p) = and V f{p) = 0. Since the decomposition 
of f(x, y) has no multiple factors, we deduce that p is a singular point of the curve 
f(x,y) = and, hence, it is a singular point of the system (JTJ). 

Our main result, Theorem |2l can only be applied when the periodic orbit 7 is 
given in an implicit way, that is, when there exists an invariant curve f(x,y) = 
such that 7 C {(x,y) \ f(x,y) = 0}. For instance, let us consider the following C 1 
system defined in all M 2 : 

x 2 ~\~ y 2 

x = (x + y) cos(x) — y(x 2 + xy + 2y 2 ), y = (y — x)(cos(x) — y 2 ) H — sin(x), (3) 

which has y 2 — cos(x) = as invariant curve. We define f(x,y) := y 2 — cos(x) and 
we have that / G C 1 (1R 2 ) and that Vf(x,y) = (sin(x), 2y). Therefore, there is no 
pel 2 such that both f(p) = and Vf(p) = 0. Moreover, f(x,y) = satisfies 
equation (0) with k(x, y) = 2y(x — y) — (x + y) sin(x). The divergence of this system 
is div(a;, y) = —^y 2 + 2 cos(x) — x sin(x) and V(x, y) = (x 2 + y 2 )f(x, y) is an inverse 
integrating factor. We denote by j n , n G Z, the oval of f(x,y) = belonging to 
the strip —ir/2 + 2irn < x < n/2 + 27m. Each oval j n , with n£Z, gives a periodic 
orbit of Q with minimal period T n > 0. The oval 70 is a hyperbolic stable limit 
cycle for system (jSJ), which can be shown just applying Theorem ^ We have, after 
some easy computations, that 



which is zero when n 7^ and it is — Air for 70. Each one of the other ovals of 
f(x,y) = 0, 7 n with ti ^ 0, belongs to the period annulus of a center as it can be 
shown from the fact that the function H(x, y) = f(x, y)(x 2 + y 2 ) exp{2 arctan(y/x)} 




—4 arctan 
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is a first integral for system J3J). Our result can be applied for any of the periodic 
orbit 7 n of this example. 

When considering a polynomial system, as far as the authors know, only alge- 
braic limit cycles are known in this implicit way A limit cycle is said to be algebraic 
if its points belong to an invariant algebraic curve. 

The paper is organized as follows. Section |2] contains the statement and proof of 
the main result of this work, i.e. TheoremEl Using TheoremEl in Sectional we will 
show that all the known algebraic limit cycles of a quadratic system are hyperbolic. 

2 Main result 

Theorem 2 Let us consider a system $IJ and 7(i) a periodic orbit of period T > 0. 
Assume that f : U C IR 2 — > R is an invariant curve with 7 C {(x,y) \ f(x,y) = 0} 
and let k(x, y) be the C l function given in (J2J). We assume that if p GW is such that 
f(p) = and V/(p) = 0, then p is a singular point of system (Q). Then, 



In order to prove TheoremEl we need to recall the definition and some properties 
of the Poincare map. Let us consider 7 a periodic orbit with minimal period T > 
for system (JTJ) and p G 7. Let C U be a neighborhood of 7 not containing any 
singular point and S = {q G W 7 | (q — p ) ■ F(p ) = 0}, where • denotes the scalar 
product between the vectors q — Po and F(p ). 

As stated and proved in pages 210 and 211 in [TSj, we have that there exists 
a 5 > and a unique function r : E — > M, which is defined continuously and 
differentiable for any q G S fl Bs(Po) such that r(p ) = and $ T ( 9 )(g) G S. As 
usual, Bs(Po) is the ball of center p and radius 5. Then, for any q G £ fl Bs(po), 
the function = $ T ( g )(g) is called the Poincare map for 7 at pq. It is clear that 
fixed points of the Poincare map, V(q) = q, give rise to periodic orbits for system 
(P). Moreover, it can be shown that V : S — » S is a C 1 diffeomorphism. 

We consider E as a subset of U C M 2 , so P is considered as a planar function 
from S C I 2 to I 2 . Hence, we notice that the derivative of ? at p , which is a 
point in S, can be represented by a 2 x 2 matrix, which we denote by DV(po). The 
following theorem, stated and proved in [2] page 118, is very useful to establish the 
stability of 7. 

Theorem 3 Let v be a non-null vector normal to F(po). Then, 







5 



Theorem |3] is proved by using the variational equations of first order related to 
system (JTJ). If $ t (x, y) is the flow related to the vector field F(x, y), we have that 
f t (D$ t (x,y)) = DF(& t (x,y)) ■ D& t (x,y) with the initial condition D& t (x, y)\ t=0 = 
I, where D means the differential with respect to the point (x, y) and / is the identity 
matrix. These equations with respect to the matrix DQ t (x,y) are the variational 
equations of first order. Since V(q) = $ r ( g )(g), the solution of the variational 
equations of first order allows the computation of DV(po) in a point po G 7. 

In order to show that the stability of 7 is determined by the value of v • DV(po), 
as stated in Theorem we consider the displacement function and we follow the 
reasoning of page 213 in [T§|. For any q G £ fl Bs(po), we have that q = p Q + sv, 
with s G (— <5/| v I , <V|v|). Since V(q) G E, we have that given s G (— 5/|v|, <5/|v|), 
there exists a <r(s) G R such that V(po + sv) = p + cr(s)v. Therefore, we have 
defined a C 1 function a : (— <5/|v|,<5/|v|) — > R and the displacement function is 
given by d : (— 5/|v(, 5/|v|) R with d(s) = - s. It is clear that d(0) = 0, 
d'(s) = cr'(s) — 1 and v ■ DV(po + sv) = a / (s)v. Since d(s) is C 1 , we have that 
the sign of d'(s) coincides with the sign of d'(0) for |s| sufficiently small as long as 
d'(0) 7^ 0. By the mean value theorem, we have that given |s| sufficiently small 
there exists a £ G (0, s) such that d(s) = d'(£)s. Therefore, if d'(0) > 0, we have 
that d(s) > for s > and d(s) < for s < 0, which implies that the periodic 
orbit 7 is an unstable limit cycle. Similar reasonings show that if cr'(0) > 1 then 7 
is an unstable limit cycle and if cr'(0) < 1 then 7 is a stable limit cycle. Theorem^ 
clearly follows from Theorem El and the fact that <t'(0)v = v ■ DV(po). 

Lemma 4 Let us consider a system Q and let f : U C M 2 — > R be a non-null 
C 1 ^) -function. There exists a C 1 function k(x, y) such that V/(g) -F(g) = k(q)f(q) 
for any q G U if, and only if, for any q EU and any t G M such that $t(g) G U, the 
following identity is satisfied: 

/(**(</)) = /(g) exp QT k($ s (q)) ds^j . (6) 



Proof. Assume that V/(g) • F(g) = k(q)f(q) for any g 6 W. We fix a point 
q E U and we define <^(t) = /($ t (g)) for any t G R such that $t(g) G U. We 
have that t belongs to an open interval (— e q , e q ) with e q > (and it may be that 
e q = +00). We have, using some of the properties of the flow and the fact /($ t (g)) = 
fc($ t (g))/($ t (g)),that: 

#f) = V/($ t (?)) " ^(?) = V/(<& t (g)) • P($*(?)) = /(**(?)) = 
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We deduce that = k($t(q)) <p(t) and <£>(0) = /(g). Solving this linear equation 

in the function cp(t) we get (pit) = /(g) exp (j^k(^ s (q))ds\. As we can consider 

the same reasoning for any g 6 W, we obtain identity (jUJ). The reciprocal is proved 
by the same reasoning. | 

Lemma 5 Let us consider a system (Jl} and 7 (t) a periodic orbit of period T > 0. 
Assume that f :U C I 2 -> 1 is an invariant curve with 7 C {(x, y) \ f(x,y) = 0} 
and let k(x, y) be the C 1 function given in (J2J). Take any p in 7. Then, 

V/(po) • I>P(po) = exp ( jf fc( 7 (*)) V/(po). (7) 



Proof. We consider the Poincare map defined in an interval of the straight line 
E containing p , V(q) = $ T ( g )(g). Since f(x,y) = is an invariant curve defined in 
U C R 2 , it is clear that for any q eW and any t el such that $t(g) G U, identity 
© is satisfied as proved in Lemma EJ Hence, 

J k{$ a (q))ds 
and differentiating this identity with respect to q we get 



/ rr(q) \ 

Vf(V(q)) ■ DV(q) = exp f J k($ s (q)) ds 1 Vf(q) 



+ 



/(g) exp I / k{$ s {q))ds 



r(q) 

(Vfc) ($ s (g)) ■ D$ s (g) ds + k(V(q))Vr(q) 



where DV(q) and DQ s (q) stand for the Jacobian matrix with respect to g of the 
functions V and $ s , respectively, in the point g. 

We evaluate the previous identity in g = p , taking into account that /(po) = 
and r(po) = an d we get identity (JJJ). | 

Proof of Theorem^ The vector V/(po) is a non-null vector that is normal to 
the vector F(p ) since /(x, y) = is an invariant curve that contains 7, and po G 7. 
The fact of V/(po) to be a non-null vector is ensured by the assumption that if 
p G U is such that /(p) = and V/(p) = 0, then p is a singular point of system 
(JTJ). Since po belongs to a periodic orbit, it cannot be a singular point. 
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Therefore, the vector v in the identity (j3J) of Theorem |3] can be replaced by 
V/(po)- Using the identity (J7J) of Lemma El we deduce that 

ex P (Y div(7(t))c^ = exp (J k{j(t))dt\, 
from which (HI) follows. ■ 



3 Hyperbolicity of the known algebraic limit cy- 
cles of quadratic systems 

We consider the families of quadratic systems with algebraic limit cycles known at 
the time of composition of this paper. These families sweep all the algebraic limit 
cycles defined by polynomials of degrees 2 and 4 for a quadratic system, as it is 
proved in j5] . In [T2"l IT3*1 IT1] , it is shown that there are no algebraic limit cycles of 
degree 3 for a quadratic system. See [S] for a short proof. In ^T], two examples of 
quadratic systems with an algebraic limit cycle of degree 5 and 6 are described. We 
study the hyperbolicity of all these limit cycles. 

The following result is due to Ch'in Yuan-shim pH] and characterizes the alge- 
braic limit cycles of degree 2 for a quadratic system. 

Theorem 6 [II Uj // a quadratic system has an algebraic limit cycle of degree 2, 
then after an affine change of variables, the limit cycle becomes the circle 

Y:= x 2 + y 2 -l =0. (8) 

Moreover, T is the unique limit cycle of the quadratic system which can be written 
in the form 

x = -y (ax + by + c) - (x 2 + y 2 - 1), , g . 
y = x (ax + by + c), 

with a ^ 0, c 2 + 4(6 + 1) > and c 2 > a 2 + b 2 . 

We summarize the four families of algebraic limit cycles of degree 4 for quadratic 
systems in the following result, which is stated and proved in [5]. We remark that 
these families were encountered previously to the work jU], but in this work it was 
shown that there are no other algebraic limit cycle of degree 4 for a quadratic 
system. System (fTUj) was first described in [T7j, system (fT^jl in [TH], system (fT4^) in 
[H] and system (fTp^) in 0. 

Theorem 7 jH] After an affine change of variables the only quadratic systems 
having an algebraic limit cycle of degree 4 are 
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(a) Yablonskii's system 

x = —4abcx — (a + b)y + 3(a + b)cx 2 + 4xy, 
y = (a + b)abx - Aabcy + (Aabc 2 — § (a + 6) 2 + Aab)x 2 
+8(a + &)cxy + 8y 2 , 

with abc ^ 0, a b, ab > and Ac 2 (a - b) 2 + (3a - b)(a - 36) < 0. 
This system has the invariant algebraic curve 

(y + cx 2 ) 2 + x 2 (x — a){x — b) = 0, 

whose oval is a limit cycle for system (jTUjl . 

(b) Filipstov's system 

x = 6(1 + a)x + 2 y - 6(2 + a) x 2 + I2xy, 

y = 15(1 + a )y + 3a(l + a)x 2 - 2(9 + ha)xy + IQy 2 , 

with < a < j^. This system has the invariant algebraic curve 

3(1 + a){ax 2 + y) 2 + 2y 2 (2y - 3(1 + a)x) = 0, 

whose oval is a limit cycle for system (|12j). 

(c) Chavarriga's system 

x = 5x + 6x 2 + 4(1 + a) a; y + ay 2 , 
y = x + 2y + Axy + (2 + 3a)y 2 , 

with ^ n+ ^2^ < a < /ias the invariant algebraic curve 

x 2 + x 3 + x 2 y + 2 a x y 2 + 2 a x y 3 + a 2 y 4 = 0, 
whose oval is a limit cycle for system ()14j) . 

(d) Chavarriga, Llibre and Sorolla's system 

x = 2(1 + 2x — 2ax 2 + Qxy), 

y = 8 — 3a — 14a x — 2axy — 8y 2 , 

with < a < -r has the invariant algebraic curve 

- + x — x 2 + ax 3 + xy + x 2 y 2 = 0, 
whose oval is a limit cycle for system (|16|). 
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In a work due to C. Christopher, J. Llibre and G. Swirszcz [II] two families of 
quadratic systems with an algebraic limit cycle of degrees five and six, respectively, 
are given. These two families are constructed by means of a birrational transforma- 
tion of system (|16p. A birrational transformation is a rational change of variables 
such that its inverse is also rational. Moreover, they prove that there is also a birra- 
tional transformation which converts Yablonskii's system (llUj) into the system with 
a limit cycle of degree 2, that is, system (0). 

The fact of the limit cycle of degree 2 being hyperbolic is stated in jTHj (see pages 
256-258) following the proof of jTUj- As a consequence, and taking into account the 
forthcoming Lemma one of the limit cycles of degree 4 (the one due to Yablonskii) 
is also hyperbolic, because this limit cycle of degree 4 is birrationally equivalent to 
the one of degree 2, as it is shown in [IT]. Our contribution is the proof of the 
hyperbolicity of the other known limit cycles of quadratic systems. 

Lemma 8 Let us consider a differential system (JI]) and a change of variables x = 
F(u,v) and y = G(u,v), where F,G are C 2 functions in U. We denote by u = 
R(u,v), v = S(u,v) the transformed differential system. Let 

dF dG OF dG 

J(u,v) :=_(„,„)_(„,„) - _(„,„)_(„,„), 

be the jacobian of the transformation. Then, 

^(F(u,v),G(u,v)) + ^-(F(u,v),G(u,v)) = ^(u,v) + ^-(u, v) + 
ox oy ou ov 

1 / dJ . . . dJ , , 

+ — r — (u, V )R{u, V) + — (u,v)S{u,v / 

J(u, V) \ou ov 



Lemma |H1 is a computational result whose proof is clear after some easy manip- 
ulations. We use it to prove the following result which states that the value of the 
integral of the divergence on the limit cycle does not change under transformations 
of dependent variables. 

Lemma 9 Let us consider a differential system (QP with a periodic orbit 7 of period 
T > and a change of variables x = F(u, v) and y = G(u, v) which is well-defined 
in a neighborhood ofj. We denote by u = R(u,v), v = S(u,v) the transformed 
differential system and by i? the corresponding periodic orbit. Then, 

f T fdP dQ\ . , xx J f T fdR dS\ /n/ x , , 
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Proof. Using the same notation as in Lemma |Hl we have that the integral 
So (faT ) ^ becomes, under the transformation of dependent variables 

x = F(u, v) and y = G(u, v), 

£ (j^(F(u, v),G(u, v)) + ^{F{u, v), G(u, v))j (0(t)) dt 

which, by Lemma |Hl equals to 

f T fdR, s dS . A , n/ „ , 

+ jf ^(w, «) + u)5(u, vfj (#(t)) dt. 

We notice that the integrand of the second integral in the former expression can be 
rewritten as d(J(u,v))/ J(u,v) and, since the change of variables is well defined in 
a neighborhood of 7, we have that this expression is a well defined, exact 1-form 
which is integrated over the closed curve so <f„d(J(u,v))/J(u,v) = 0. | 

Therefore, in order to prove that all these families of limit cycles are hyperbolic, 
we only need to study the stability of the limit cycles of systems ljI2|L (jH|) and 
The hyperbolicity of the two limit cycles described in ^T] is shown by the fact that 
they are birrationally equivalent to (|16|). 

Theorem 10 Eac/i one of the limit cycles of systems Q12j). (|I%J) and flTtjj) hyper- 
bolic. 



Proof. In order to prove the hyperbolicity of the limit cycles of systems ()12|). 
()14|) and (JTOJ) we use the same process for all of them. These systems depend on a 
parameter a which belongs to a certain open interval when the limit cycle 7 exists. 
We denote by T > the period of the limit cycle and by T>(a) the value of the 
integral div(7(i)) dt for any value of the parameter for which the limit cycle 
exists. This value decides the hyperbolicity character of the limit cycle 7 in the 
system with parameter a. By virtue of Lemma El we may consider any birrational 
transformation of these systems well defined in a neighborhood of the limit cycle 
and we may consider the transformed system instead of the previous one because 
the value of the integral f Q div(7(t)) dt does not change. 

Using Theorem we have that: 

V(a)=f div(7(t))dt + w ( [ div(j(t))dt- ( k(j(t))dt 
Jo \Jo Jo 
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where k is the cofactor of the invariant algebraic curve containing the limit cycle 
and w is any real number. 

We show that the function T>(a) has no zero when a belongs to the interval of 
existence of limit cycle by choosing an adequate idgR and parameterizing the limit 
cycle 7. The way of choosing the adequate value of w is purely heuristic, although 
we expect that this choice is related to some geometric property. We find it very 
surprising that it is possible to choose w = — 3 for each one of the three families of 
systems. 

Hyperbolicity of the limit cycle given by the algebraic curve ()17|) for system (fTEJ). 

The stability of the limit cycle 7 is given by the following function of the pa- 
rameter a of the system, T>(a) := f Q div(7(t)) dt, where div(x, y) = 2(2 — 5ax — 2y) 
is the divergence of system (j!6|) and T > the period of the limit cycle. Theorem 
12] gives: 

div(7(t))d*= / k(r/(t))dt, 
'0 Jo 

where k(x, y) = 4(2 — 3ax + 2y) is the cofactor of the invariant algebraic curve (|17|). 
So, given any real number w, we have that: 

V(a) = [ div(7(t))dt + w [ (div - k){j(t)) dt 
Jo Jo 

r-T 

((1 + w)div - wk) (j{t))dt. 



We consider the following parameterization of the oval of the algebraic curve ([17)1: 



X{T)=T , fc(r) = iiMS Eg, (19) 

It 



where T\ = 1 ^ 4tt , r 2 = 1+ ^ 4a and the parameter r 6 (ti,^). The positive 
sign y+(r) gives a half of the oval and the negative sign y~{r) the other half. One 
of the endpoints of both parameterizations is (x\,y\) = ( 1- ^~ 4a , — 1+ Vi-4 a^ anc j 
the other endpoint is (x 2 ,y 2 ) = ( 1+ ^P° , ^i±^0«). We have that the vector field 
in (xi,yx) is (0, 6\/l — 4a) and in (#2,2/2) i s (0, — 6\A — 4a), so the flow on the 
limit cycle is clockwise. The line 2ax = 1 cuts the limit cycle in two points with 



ordinates ±y ~ a -, which are given respectively by y±(l/2a). We have the 
following relation between the differentials: dr = P(x(r), 2/±(t)) dt where P(x, y) = 
2(1 + 2a; - 2ax 2 + 6xy). Then, 

X>(a) = / ((l + w)div- wk) (i(t))dt 
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T2 



((1 + w) div — wk) 
P 



( r > y+( T )) dr + 



((1 + w) div — wk) 



P 



(r, y_(r))dr 



fT2 


[( 


1 n 





((1 + w) div — wk) 
P 

((1 + w) div — wk) 
P 



( r > y+( r )) 



dr. 



For w = — 3 and substituting by the parameterization, we get, 

V / ar(r^Ti)(T^ : r) 



Via) = 8 



r(l + 8r + ar 2 



■dr. 



Since r x > and r 2 > r x for any a G (0, 1/4) and the integrand v r( - 1+8r 



•y/ ar(r-ri)(T^2— r) . 



(l+8r+ar 2 ) 



IS 



also strictly positive and well defined for any r G (ri, r 2 ) and a G (0, 1/4), we have 
that P(a) > for all a G (0, 1/4), which implies that the limit cycle in system (fT7)|) 
is hyperbolic (and unstable). 

Hyperbolicity of the limit cycle given by the algebraic curve ()13j) for system (|12p. 

In order to simplify our computations, we consider the following birrational 
change of the parameter, a = 3c/ (4 + 5c), with inverse c = 4a/ (3 — 5a) and we have 
that c G (0,1/2). 

We consider the birrational change of variables (x, y) — > (u, v) given by (x, y) = 
(X(u,v),Y(u,v)), where 



X(u,v) 
Y(u,v) 



2(1 + 2c) 
c 2 (1 + m) 3 
12(1 + 2c) 2 



c - 2(1 + c)u + cu 2 - 2Vl + 2cv), 

c — 2(1 + c)u + cu — 2 VI + 2cu). 



c 2 (4 + 5c)(1 + m) 4 

The inverse of this change is given by («, f ) = (U(x, y), V(x, y)) with 

v = 6(l + 2c) X _ 
(4 + 5c)y 



V 



VT+2c~ 



1 + 2c)(54cV + 18c(4 + 5c)x 2 y - 6(4 + 5c) 2 xy 2> 



(4 + 5c) V 

The jacobian of this change of variables is: 

dU dV dUdV 324 c 2 ( 1 + 2c) 5 / 2 x 4 
dx dy dy dx (4 + 5c) 4 y 5 
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+ vTT2c. 



which can be seen to be well defined and different from zero in all the points of the 
oval of the curve given in ()13j). 

We get a transformed system in which we reparameterize its time t multiplying 
by c 2 (4 + 5c) (1 + u) 3 /(12(l + 2c)). This reparameterization does not affect the 
direction of the flow on the limit cycle. The new system reads for: 



u = -2u (cu + 4 + 9c) (cu 2 -u + c) - 2y/l + 2c (cu 2 - (4 + 5c)u + 2c)v, 
v = -c 2 v / TT2^(m + 1) 2 (m- 1)(3u + 2) (20) 
-(cm + 4 + 9c)(3cw 2 - 2u + c)v + 2^1 + 2c (4 + 5c - 3cu)v 2 , 

and the limit cycle is transformed to the real oval of the curve v 2 + u(cu 2 — u + 
c) = 0. This algebraic curve is invariant for system (|2(Jj) with cofactor k(u,v) = 
4a/1 + 2c (4 + 5c — 3cu)v — 2{cu + 9c + 4) (3cu 2 — 2u + c). The divergence of system 
Pty is div(«, i>) = 2v/IT2c" (12 + 15c - 8cu)v - [llcV 5 + c(28 + 81c)u 2 + (5c 2 - 
54c-24)u + 3c(4 + 9c)]. 

We consider the following parameterization of the oval of the algebraic curve 
v 2 + u(cu 2 — u + c) = 0: 



uir = r, 



v±{t) = ±^ct{t-t x ){t 2 -t), (21) 



where t\ = 1 ~ v 2 1 ~ 4c2 , r 2 = 1+v/ 2 1 ~ 4c2 and the parameter r G (ti,t 2 ). We notice 
that for c G (0,1/2), we have that < T\ < 1 < r 2 . The endpoints of both 
parameterizations are (t 1? 0) and (t 2 ,0), and the vector field at the point (t;,0) is 
(0,cVl + 2c(Ti + l) 2 (3ri + 2)(1 - n)), for i = 1,2. Therefore, the flow on the 
limit cycle is clockwise. We follow analogous arguments to the previous example to 
deduce that: 



V(c) 



((1 + w ) div — wk) 
P 

((1 + w) div — wk) 
P 



( r ? v+(t)) 



(r, u_(r)) 



where P(u,v) is the polynomial which defines u = P(u,v). For u; = —3 and 
substituting by the parameterization, we get 

v / Z" 7 " 2 \/ct(t — Ti)(t 2 — r) 

X> c = Sv/TT^ / t N v , , — r dr. 

KJ J T1 (r + l)(cr 2 + (17c + 8)r + 4 + 8c) 

Since < T\ < 1 < r 2 for any c G (0, 1/2) and the integrand is strictly positive 
and well defined for any r G {ti,t 2 ) and c G (0, 1/2), we have that T>(c) > for 
all c G (0, 1/2), which implies that the limit cycle in system (|2*0|) given by the real 
oval of v 2 + u(cu 2 — u + c) = is hyperbolic (and unstable). Hence, using Lemma 
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we have that the limit cycle given by the oval of the curve (|13|) in system (|12|) is 
hyperbolic and unstable. 



Hyperbolicity of the limit cycle given by the algebraic curve (frSj) for system (fPfj) . 

We consider the following birrational change of variables (x, y) — > (u, v) with 

(x,y) = (X(u,v),Y(u,v)) where 

(X(u,v),Y(*,v)) (_2 " 2,,! 



v + 1 + u + 2av? 
The inverse of this change is (u, v) = (U(x,y),V(x,y)) with 

TTf \ y \r/ \ y 2 (y + 2) 

x x z x 

This change of variables is well defined in a neighborhood of the real oval of the 
curve (fT5|l and its jacobian is 

dUdV _dUdV 2_ 
dx dy dy dx x 3 

The algebraic curve (fH)j) is transformed to v 2 +4au 2 (u— 1) — (u+1) 2 = 0. We consider 
the transformed system in which we make a reparameterization of its time t which 
consists on multiplying by v + 1 + u + 2au 2 . This reparameterization reverses the 
direction of the flow on the transformed limit cycle and the new system is written: 

u = (u + l) 2 - 4au 2 (u - 1) + (1 - 3u)v, 

(22) 

v = 2(u+l)(3 + u + 2au-au 2 ) + (l + Aau + u-Qau 2 )v-5v 2 . 

The divergence of this system is div(u,v) = 3(1 + u) + 12cra — 18cra 2 — 13t>. The 
algebraic curve v 2 + Aau 2 (u — 1) — (u + l) 2 = is invariant for system (}2*2*|) with 
cofactor k(u,v) = 2(1 + u + Aau — 6au 2 — 5v). The real oval of this curve is a 
hyperbolic limit cycle for system (j22j) if, and only if, the real oval of the curve 
()15|) is a hyperbolic limit cycle for system (JHJ), by Lemma EH We are going to 
compute the function V(a) = / Q T div(7(t)) dt, where 7 is the real oval of the curve 
v 2 + 4au 2 (u — l) — (u+ l) 2 = 0. To do so, we parameterize this oval by: 

U ( T ) = r , v± (t) = ±y/Aar 2 {\ - r) + (r + l) 2 , 

where r takes values between r± and r 2 . The values T\ and r 2 are the two small- 
est roots of the polynomial g(a,r) := 4ar 2 (l — r) + (r + l) 2 in r. We con- 
sider a in the interval between (17\/l7 — 71)/32 and 0, which are the values of 
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the parameter for which the limit cycle exists. Since the coefficient of the high- 
est order term of g(a,r) is —4a which is strictly negative, g(a, — (3 + \/Tf)/2) = 
4(29 + 7VTT) [a - (17VV7- 71)/32] > and o(a, -1) = 8a < 0, we deduce that 
7~i < —(3 + \/Tf)/2 < r 2 < —1. We denote by P{u,v) the polynomial which defines 
m = P(u,v) in system (|2*2*|) . We consider the point with coordinates: 



(U ,V ) 



3 + V17 




and we have that Vq — g(a,uo) = and P(wo, t>o) = vq[vo + (11 + 3vl7)/2]. We 
deduce that P(u , v ) > and that the flow on the limit cycle is clockwise. Using 
analogous arguments as in the previous examples we conclude that: 



V(a) 



n 



((1 + w) div — wk) 
P(u,v) 
((1 + w) div — wk) 



T. 



,v + (r)) 



(r, u_(r)) 



P(u,v) 

For w = —3 and substituting by the parameterization, we get 

\/g{a,r) 



dr. 



V a] 



(r-l)r(ar + 2) 



We have that n < -(3 + Vl7)/2 < r 2 < -1 and that r - 1 < and r < < -2/a 
for any r G (ti,t 2 ) and for any a G ((17VT7 — 71)/32,0). Hence, the integrand is 
strictly positive and well defined for any r G (ti,t 2 ) and a G ((17\/17— 71)/32,0). 
We deduce that V(a) > for all a G ((17>/l7 - 71)/32, 0), which implies that the 
limit cycle in system (|22j) given by the real oval of v 2 — g(a, u) = is hyperbolic 
(and unstable). Thus, using Lemma El and the sign in the change of time, we have 
that the limit cycle given by the oval of the curve (J15)) in system (114)) is hyperbolic 
and stable. | 



Appendix 

The aim of this appendix is to present some relations among elliptic integrals which 
the authors obtained by using the identity given by Theorem El for systems (|T2*|) and 

(Hi). 

Before the presented proof of Theorem E3 the authors got its proof for systems 
(JT2J) and (fT7)|) by computing the corresponding integrals which give place to elliptic 
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integrals. The identity given in Theorem |2] was used to encounter a Fuchs equation 
for the function T>(a). After some thorough analysis of this Fuchs equation, we 
deduce the non- vanishing of the function T>(a) for any value of the parameter in 
which the limit cycle exists. We are not going to give this proof, but we think 
that the relations among elliptic integrals obtained by the former reasoning are 
interesting by themselves. Hence, we give the identities obtained which, as far as 
we know, do not appear in any book of tables of integrals and relations between 
classical functions. On the other hand, we also give the obtention of the Fuchs 
equation for the function T>(a) in the case of system (jTBj) . 



Identities among elliptic integrals 

The functions involved in this subsection are the complete elliptic integrals of first, 
second and third kinds, denoted by K(o>), E(a>) and H(k, uj), respectively. We recall 
the definition of these functions: 



K(«) 



n/2 de r 1 dt 



l-ujsm 2 (9) Jo -t 2 )(l-uJt 2 y 



-tt/2 i pi Pi _ , 2 

E(w) = / \ l-ujsm 2 {6)d6 = / V . dt 



h(k, to) 



v/2 de r 1 dt 



(l-Ksin 2 (0))v/l-a;sin 2 (0) Jo (1 - Kt 2 )y/(1 - t 2 )(l - ut 2 ) ' 
and their derivatives: 

k'H = 9n 1 , em-Ikm, 

2(1 — uj)uj zu 
E'H = i.E( W )-^K( W ), 

dU(K, Uj) 1 w . 1 . K 2 — U T-r/ x 

ok 2k[k — 1) 2[k — — K) 2k{k — l)(o; — k) 

uj) 1 s 1 . - 

y = T( w 7T E(u;) + — rn «, W 

OOJ 2{k — uj){uj — 1) 2{k — uj) 



We use the following parameterization of the oval of the algebraic curve (|17j) to 
explicitly compute the integrals for the system (|16j). We parameterize the oval by: 

1±V [-Mr-^ )Ml i 
It 
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where T\ = 1 ~^~ 4a , r 2 = 1+v ^~ 4a and the parameter r G (ti,t 2 ). The positive 
sign y+{r) gives a half of the oval and the negative sign y-{r) the other half. The 
explicit computation of the integrals for the system (|T6*j) gives that the identity Q 
stated in Theorem |2] reads for: 

- 9 K(w ) + c + n(cj + , cj ) + c_ n(cj_, cj ) = 0, (24) 

which is valid for a G (0, 1/4), where 

2y/l - 4a 2V1 - 4a 9 - Vl - 4a 



tun = , cj+ = . , c+ = ± yl6 — a. 

1 + VI -4a 9 + VI -4a±2V16 - a 2 



The derivative of the expression in (j24j) with respect to a gives place to the same 
identity (|24p. In fact, when computing the derivative with respect to a of the 
expression given in ([24)1 . using the described formulas of derivation for these elliptic 
integrals, we get — 1/(1 — 4a+ Vl — 4a) times the same expression (J2*3|) . This simple 
factor is different from zero when a 6 (0, 1/4). 

In the same way, we can explicitly compute the integrals involved in the identity 
(J3J) stated in Theorem El for the system (|12|). via using the parameterization of the 
oval of (fT3j) given by: 

2(1 + 2c) 



c (1 + t) 

y ± (r) = c 2^t)(l C +T)* ic " 2T " 2CT + Cr ' ± 2v^raV(-c + r-c7-2)), 

(25) 

where r G (n, r 2 ) with n = W 2 ~ 4c2 and r 2 = 1+v 2 1 ~ 4c " . It is clear that < n < 
1 < r 2 for c G (0, 1/2). We define g{c,r) = —t(c — r + cr 2 ) = cr(r — ti)(t 2 — r), 
which is strictly positive for all r G (n, r 2 ). The explicit computation of the integrals 
involved in the identity (@J) gives: 

5 Kfco) + n(?+, ? ) + CL n(?_, ? ) = 0, (26) 
which is valid for c G (0, 1/2), where 

2V1 - 4c 2 2Vl - 4c 2 

?o = , „ , ?± 



1 + Vl - 4c 2 ' =c 9 + 17c + Vl - 4c 2 ± V64(l + 2c) 2 + 
-2(24 + 47c ± 3 V64(l + 2c) 2 + c 2 



c 2 ' 



9 + 17c + Vl - 4c 2 ± V64(l + 2c) 2 + c 2 

The derivative of the expression with respect to c gives place to the same 
identity flggp. 

The authors have not been able to give an analogous identity related to system 
(fTlJl due to the fact that the corresponding integrals require much more computa- 
tions to be identified with the elliptic integrals. 
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Fuchs equation for T>{a) in system ( 1161 ) 



In this part of the appendix we develop the way we obtained a Fuchs equation 
for the function T>(a) in system ([16)1. via using the relation (j2U). We think that 
the fact of obtaining a Fuchs equation satisfied by this function is interesting to 
further understand the stability of algebraic limit cycles for polynomial systems. We 
obtained a similar Fuchs equation for system (JT2J), but we do not state it because the 
equation itself does not give any further information about the properties of system 
(H2J) and the way it was obtained is completely analogous to the way equation (J2*T|) 
for system (fTBj) is obtained. 

Let us consider system (J16j) and we parameterize the oval which contains the 
limit cycle by (|23|). Taking the notation described in the previous subsection: 



UJQ 



and 



2V1 -4a 



1 + VI -4a 



2\A - 4a 



9 + VI - 4a ± 2y/W 



c± = ± vlo — a, 



/i = ^1 + Vl-4a, b± = 2(4 ± VlQ -a) c±, 
we explicitly compute the value of V(a): 



V(a) 



div(7(t))dt 



V2 



fi\/16 — a 

We compute the successive derivatives of £>(a): 



[-34V16 - aK(cj ) + b + Il(u + , oo ) - 6- n(w_, w„)] 



D'(a) 



-4^2 



, ...... .. . Vl6-aK(o;o) + 

\i (16 — a)^ L 

-C + L T (co' + ,Co'o) + C_n(tU_,Co'o) 



2/x 2 Vl6-a 



E(w ) + 



6^ 



/x(16-a) 5 /2 



(10a 2 + 33a-64) v / 16^ 



3a(l - 4a) 



K(w ) + 



(73a 2 - 420a + 128)/i 2 v / 16" 



6a 2 (1 -4a) 



■E(a>o) + c+n(u; + ,^o) - c_il(^_,^ ) 



P'"(a) 



-V2 



/i(16 - a) 7 /2 



180a 4 + 1347a 3 - 9685a 2 + 25664a - 4096)^16 - a 
a 2 (l -4a) 2 



K(w ) 



(1812a 4 - 20259a 3 + 102164a 2 - 60544a + 8192)/i 2 v / 16 - a . . 
+ 2a 3 (l-4a) 2 ^ + 

— 15c + U.(u + , u> ) + 15c_n(u;„, u ) 
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By elimination of independent functions and using the identity (|24j) we obtain the 
following third order homogeneous differential equation of Fuchs type for T>(a): 

8(a - 16)a(4a - l)(17a + 8) 22'"(a) + 4(612a 3 - 4119a 2 - 2600a 

+512) V"{a) + 6(a - 2)(289a + 528) 22'(a) + 3(17a + 64) 22(a) = 0. ( 27 ) 

An easy computation shows that 22(1/4) = 0, 22'(l/4) = -8V2n/9 and 22" (1/4) = 
98v / 2tt/27. Hence, equation ()27)1 univocally determines the function 22(a) defined 
in a G (0, 1/4]. A thorough analysis of the properties of 22(a) gives that 22(a) > 
for a e (0,1/4). 

We remark that using identity (}2"4"j) we get a Fuchs equation of order 3 for 22(a). 
If we did not have this relation, we would get an equation of order 4, which would 
make the analysis of properties much more difficult. We notice that this Fuchs 
equation is an interesting alternative method to prove the hyperbolicity of the limit 
cycle in system JTEJ). This kind of equation may exist for all algebraic limit cycle of 
a planar polynomial system and may let distinguish its hyperbolic character. 
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